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SUMMARY 



Computat ions* were made to determine the effect of 
some airplane design trends on the fin area and the dihe- 
dral angle required for lateral stability. The specific 
factors studied were wing loading, moments of inertia in 
roll and yaw, wing chord, and tail length. The airplane 
for which the computations were made was chosen to be 
representative of pursuit-type airplanes, hut the conclu- 
sions drawn are applicable to any type of airplane char- 
acterized by the parameters and the data employed. 

The results of the computations are presented in the 
form of diagrams of variations of fin area with dihedral 
anglefor neutral stability. The results indicate that 
increasing the values of wing loading and moments of iner- 
tia makes the attainment of lateral stability increasingly 
difficult and even doubtful in some cases for airplanes 
with flaps deflected. The fin area required for lateral 
stability is more dependent on fuselage size than on wing 
area. _ . w 

INTRODUCTION 



In present airplane design, the trend is to use high- 
lift devices, high power, and weight distributed along 
the wings. Since. the result of this. trend is to increase 
the wing loading, the radii of gyration in .roll and yaw, 
and the fuselage size with respect to the wing area, the 
amount of fin area and dihedral angle required for lateral 
stability will undoubtedly be affected. 

The amount of fin area and dihedral angle required 
for the lateral stability of an airplane may be computed 
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from the equat ions:. o~f njotion ; fO;r *sm.a-ll oscillations, pro- 
vided the aerodynamic characteristics of the airplane are 
known. Although few aerodynamic data are available for 
some of the, parameters involvedv^seful information can 
"be obtained by assuming aerodynamic characteristics and 
studying the resulting -trends. Some investigations of 
this nature have been reported in references 1 and 2. 
The data given in both references are limited in scope, 
and the direct effect of changes in the fin area and the 
dihedral angle is difficult to visualize. Also, the data 
given in reference 2 do npt include the conditions for 
high wing loadings and large radii of gyration. 

The present report covers' the results of an investi- 
gation planned to show the. general effect of changes in 
wing loading, radii of gyration, wing chord, and ta'il 
length on the amounts of f in : area and dihedral angle re- 
quired for lateral stability. The changes in wing load- 
ing and radii of gyration are intended to represent the 
present trends in pursuit-type airplanes. The. results, 
however, apply to any. type of airplane characterized by 
the parameters and the aerodynamic data used. 

Computations were made for neutral spiral and oscil- 
latory divergence, using representative values of the 
stability derivatives for the airplane as well as for 
changes in fin area, dihedral angle, and aspect ratio. 
The results of , the computations are given as diagrams of 
dihedral angle against fin area for neutral spiral and 
oscillatory divergence. 

The airplane parameters used in this investigation 
were the wing loading, radii of gyration, wing chord, and 
tail length. The variations in wing loading and radii of 
gyration were planned to represent present-day trends in 
pursuit-type airplanes. One parameter was varied at a 
time while the others were kept at a mean value. In some 
cases a change -in one parameter is either directly or in- 
directly the result of a. change in other .parameters. The 
following table gives, in each group, all. changes caused 
by .the change in a p^rt icuiar -parameter . 
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where 



W/ S 



kx/b and k z /b 



wing loading ■ 
weight of airplane 
wing area 

ratio of airplane density to air den- 
sity computed for wing span of 40 
feet and standard air density at 
sea level (m/pS w b) 

ratio of radius of gyration about re- 
spective airplane axes to wing span 
(These valxies have been spotted on 
a plot (fig.l) of values of k x /b 
against k z /b from available design 
data of pursuit-type airplanes,) 

ratio of distance from center of grav- 
ity of airplane to center of pressure 
of fin to wing span 
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A. 


aspect ratio of wing 


m 


airplane mass 


P 


air .density- 


S f 


fin area (total vertical tail area) 




wing span 



STABILITY DERIVATIVES 



Because few data are available for some of the aero- 
dynamic characteristics required, particularly for the 
interference between component parts of the airplane, two 
basic assumptions are necessary for computing the lateral- 
stability derivatives. These basic assumptions are: 

1. The stability derivatives for each component part 
of the airplane are added to give the value for the com- 
plete airplane, that is, the interference effects are asr 
sumed to be zero. In some cases the interference effects 
are large and depend upon the arrangement of wing, fuse- 
lage, and tail-surface combinations. (See references 3 
and 4 for some of the interference effects caused by wing 
location on the variation of the derivatives that depend 
upon sideslip.) 

2. Neither the propeller nor the slipstream has any 
effect on the stability derivatives. These effects may 
be large for some of the derivatives. 

The stability derivatives used and their variations 
with. the airplane parameters were computed by the follow- 
ing relations, which are intended to include the effects 
of the parameters and particularly the changes in fin area 
and dihedral angle. The stability derivatives are the in- 
stantaneous rates of change of the aerodynamic coeffi- 
cients with attitude or angular velocity when the attitude 
angle or the angular velocity is zero. For convenience 
the derivatives are written in the form c l ft > c l » ' $\ > 
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.. .dCi do; 
etc. instead of dC^db, ~-^r. ~7rty etc - where rep- 

. d (2v) d ( 2 "v). 

resents the force or moment coefficient and p and r 
represent the angular. rates of rotation, rolling and yaw- 
ing, respectively. 

The relations for the derivatives are: 

Cy " = Cy . + C v + Cv 

y p(fuselage) Y £(wing) Y (3(fin) 

= - 0.020 A + 0 - 3.48 It ( \) 

where ^ 

Cy force coefficient along Y axis of airplane 

$ angle of sideslip, radians, positive when the 
right wing is into the wind 

and the subscripts y indicate the contribution of the cor- 
responding Tsart i0 f the airplane to C v . (Most wind- 

tunnel data have been given as dCy 1 /^' 1 where Cy 1 is 
the force along the Y wind axis and is the angle 

of yaw, degrees (57.3 p = - \|/» ) • therefore C v 

dC Y ! 

= - Tfr 57 • 3 ~ V. 

The constant' 0 ,020 in Cy = - 0.020A was 

£( fuselag.e) 

computed from data given in' reference 3. The fact that, 
in this investigation, the fuselage size remains constant 
and the wing area varies inversely, as the aspect ratio 
necessitates the use of this expression. 

Although Cy varies with dihedral ( references 

(3(wing) 

3, 5, 6, and 7), this derivative has not been included be- 
cause it is counteracted, at least partly, by the deriva- 
tive of the side force due to rolling C- . 
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The rate of change of normal force on the fin with. 

S 

angle of sideslip, C_ , is equal to - C 3.48 

Y P(fin) . V ■ S w J 

in terms of the wing area. The term Sf/S. w is the ratio 

of fin area to wing area (a variable for this investiga- 
tion) and 3.48 is the rate, of change of normal force co- 
efficient on the fin with sideslip angle p for an ef- 
fective aspect ratio of 3. The value obtained with the 
fin used on the model in reference 3 is about 3.48. 



where 

rolling- moment coefficient 

c ]o/ . = . K x r 

L P(wing) 1 - 

Z x varies with aspect ratio; it is equal to 
- 0.0175 for A = 16, - 0,0141 for A= 8, 
and - 0.0114- for A'= 4. These values 
of Ki were obtained from data given 
in reference 6 for a wing of a 2:1 taper 
ratio . 

T effective dihedral angle of wing, degrees. 

Effective dihedral angle is used through- 
out this report as a fictitious angle 
that would give the wing the value of 
C 7 . . Wing plan form and elevation, as 
P ~ : 

well as large interference effects, con- 
tribute various amounts of C^. (See 

references 3, 5, and 7.) 
The contribution of the fin to the value of is 

where ^5 - ll S i n i s the assumed ratio of the distance 
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of the center of pressure on the fin above the horizontal 
plane to the wing span* The value of x was assumed to 
"be the distance from the fuselage center line to the cen- 
ter of pressure of the fin. Upon the assumption that the 
distance from the top of the fuselage to the center of 
pressure of the fin varied with the square root of the 
fin area and with the use of data from reference 3, the 

relation | = 0.025 + 0 .23 JsJJl^ was obtained and used 

for all variations of the fin area. The value of i, 
the angle of fuselage center line to the horizontal, was 
zero for all cases when the lift coefficient C L was 

equal to 0.2. For other values of , the value of i 

depended upon both C L and A. The value of i for C L 

== 2,8 was the same as for C T = 1.4. The increase in C T 

lr L 

was obtained by merely deflecting the flaps with no change 
in the angle of attack. The term 3.48 S f /s w was used in 

the expression for C v 

Y P(fin) 

% " ° n f3(wing) + C ap(fuselage) + C np( fin ) 

= K 2 - 0.009A + ~ ^3.48 ( 3 ) 

where C n is the yawing-moment coefficient. 

A value, of. K 2 of 0.009 was used for plain wings 

and of 0.030 when the flaps were deflected. Values for 
K 2 may vary considerably with wing forms. (See refer- 
ences 5 and 7. ) 

The value C n = - 0.009A was obtained from data in 
P 

reference 3 and may be expressed as a function of A be- 
cause the fuselage and the wing span are constant with 
change of wing area. 

The contribution of the fin to the va of C i« 

the computed variation of the normal force with 0 on the 
fin C *{3(fin) times the nondimensional lever arm l^/b. 
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07=0, = K 3 (4) 

l p l p(wing) 

where the values of K 3 (from reference 6) are -0.6 for 

A = 16, -0.5 for A = 8, . and - 0.4 for A = 4. The 

values of Ci - and C7 , x are probably 

*p(fin) ^(fuselage) 

very small as compared with C lp( w i ng ) an<i therefore 
have not "been used. 



+ 2 



where K4 varies with the aspect ratio and is (from ref- 
erence 6) - 0.089 for A = 16, - 0.065 for A = 8, and 
- 0.040 for A = 4. 

The expression for C~ , v , that is, 

n p(fin) 

2 ^1 It sin i^ ^3.48 is the rate of change of 

the yawing moment with the rolling velocity, produced "by as- 
suming that the normal force on -the fin is proportional 
to the ang] e induced at. the center of pressure of the fin 
by the rolling velocity; actually the induced angle is a 
variable along the fin. 



07 =01. x + C ? \ = 0.25 C T 

fc r ^r(wing) Ir(fin) 



where the value of 0.25 was obtained from data given in. 
reference 6, and the term for C7 is of the^'same 

Mfin) 



NACA Technical Note No. 814 



9 



form as for C n , „ x ; but, in this case\the angle is 
p(fin) 

induced by the yawing velocity and the value of C7 

°r(fin) 

therefore is probably more nearly representative of the 
actual value of the derivative than the value given for 

Cn p(fin)' 



°n r - C n r ( fuse iage) + ° n r(wing) + C *r(fin) 

= - 0.00125A + (e 5 C l 2 - K 6 ) - 2 (^.) 2 ( 3 ,48 ^ (7) 



where Z5 (from reference 5) is - 0.0113 for A = 16, 
- 0.02105 for A = 8, and - 0.03838 for A = 4. The con- 
stant K 6 depends on the profile-drag coefficient of the 
wing (reference 6) and is assumed to be 0.003 for the 
plain wing and 0.030 for the wing with flaps deflected* 
The value - 0.00125A was assumed for the fuselage and 
varies inversely with the wing area. 

The expression for C n r ( fin )» that is, 

( l t V / S f \ 

2^~^/ ^3.48 ~j , is the nondimensional form of the 

rate of change of yawing moment due to the fin with yaw- 
ing velocity. 



THEORETICAL EQUATIONS 

The equations used to compute the boundaries of neu- 
tral spiral and oscillatory stability were developed from 
the theory of small oscillations, as given in reference 2, 
simplified for the level-flight condition. The equations 



are of the form; 



AA 4 + BX 3 + CA 2 + DA + E = 0 
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where 



A = 1 
E = - Y v - n r - t 



P 



and 



C = l p n r - l r n p + Y v (l p + n r ) + |0.n v 

D = Y v ^l r n p - l p n r ^ + [Ll y (n p - ^ - ^n v l p 

E = I" ( |ll v n r - ^ n v l r) 



^ ' "( VO 2 " 



2 1X011 p 



txn v = --— ---- g 



1 p = (v";y 



P / , N2 



( k z/*y 
(k x / 0 y 



. . 4 *r 



1 c 

n r = ~ 



(V»)' 
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The necessary, conditions for' stability are that the' 
values of A, B, C» D, E, and Routh*s discriminant *R . 
shall all be positive. Since the stability derivatives- 
and the airplane parameters used for this investigation 
are such that the values of A,.B, and C are positive 
and the values of D are positive, when R is positive 
for any practical value. of effective dihedral angle, sta- 
bility depends on the values of E and R, When the 
value of E becomes negative, the airplane becomes spi- 
rally unstable. The la t eral 0 sc il lat ians increase in 
amplitude when the value of R becomes negative. The 
limits of the stable region are therefore defined by the 
values of T and Sf/s w when 

E = 0 

and 

R = BCD - D 2 - B 2 E = 0 



RESULTS AND DISCUSSION 



The results are presented in the form of diagrams 
(figs. .2 to 5) showing the variations in the computed 
boundaries of spiral and oscillatory stability with 
Sf/S w . (the ratio of the effective fin area to the wing 

area) and with I\ the effective dihedral angle. Figure 
6 is a replot of data given in figure 4(c) with the ordi- 
nate given as. (s f /S w J ( 8/a) . . Figure 7 is a replot of 

the data given in figure 5(c) with the ordinate given as 
( s f/ s w) \}t/^) 9 (effective tail volume). The value of 

Sf/S w where C- " is zero, the value of C n required 

P p 

for'neutral weathercock stability, is indicated in figures 
2 to 7. ' 

The results, in general, indicate that the value of 
r required for spiral stability increases with C^ and 

S f /S w , The value, of ;Sf ; /S w required for oscillatory 

stability decreases iwith r ; f or several degrees and then 
increases with continued increase in* T. For some cases 
the diagrams do not cover a sufficient range of F to 
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show the increiase in Sf/S w . The rate of change of Sf/S w 

w.ith F required for oscillatory stability increases with 
and with some of the airplane parameters. 

The stability "boundaries for = 2.8 may appear to 

be inconsistent with the other lift coefficients. The 
differences, however, result from the changes in C^, 

Cnp» an( * Gn T produced "by the flaps and from the assump- 
tion that the angle of attack is the same as that with no 
flaps, that is, = 1.4. 

Effect of wing loading on stability boundaries. - The 
value of jju , relative density of the airplane to air, has 
no effect on the boundary of neutral spiral stability. 

(See fig. 2 and the expression E =C L /2 (\xl^n T - |xn v l t ) = 0 . ) 

The change in the boundaries with a variation in re- 
sults from the variation of the derivatives with and 
not from the presence of \x in the expression for B. 

For oscillatory stability, increasing \x increases 
the value of Sf/.S W required for a given value of T 
and decreases the. value of T for a given value of S f /s w . 

This effect is small for small values of T but becomes 
large when F is increased more than a certain amount and 
is particularly important for high values of and 

(See figs. 2(c) and 2(d),) 

The resulting effect of increasing (j,,* is. to increase 
the difficulty of obtaining lateral stability, especially 
with large values of C L and with flaps deflected. 

Effect o f ra d i i of gyration in .vol 1 _and _yaw on sta- 
bility boundaries . - The effects of k x /b and k z /b ( fi S* 
3.) are practically the same as the effects of |ju The 
values of k x / b and k z /b used have a greater effect on 

the oscillatory boundaries than have the values of 
There is no effect of k x /b and kg/b on the spiral 

boundar i es . 

The data given in figures 3(c) and 3(d) indicate 
that, for large values of C^, it is impossible to sat- 
isfy the two conditions both with and without flaps. The 
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value of Cj,p represented in this report as an effective 

T usually increases with C L if the wing has sweepback 

or certain tip shapes (references 5 and 7). Although this 
variation in effective T with C L might he useful in 

obtaining "better stability characteristics, it may cause 
serious oscillatory instability at the high values of C L 

with large values of \± or of k x /b and k z/ b * 

^l£Qc^_Q-f,_changing wing chord o n stability boundaries. 

Changing the wing chord without changing the span, the 
fuselage, or the total weight of an airplane changes A, 
(ji , and S w . Because the fin area, as used in this report, 

is represented as a ratio of the wing area, S f /s w , the 

actual fin area is inversely proportional to A. These 
factors must be considered when the effects of changes in 
wing chord, given in figure 4, are compared. Figure 6, 
which is included to show the effect of these factors, is 
a replot of figure 4(c) based on the actual value of the 
fin area. 

Decreasing the wing chord (increasing A and |a) , 
in general, decreases the value of T necessary for spi- 
ral stability (fig. 4). This effect must be due to A 
and to the ratio of wing area to fuselage size because, 
as previously explained, |i has no effect on these 
changes. When the actual values of S f are considered 

(fig. 6), the variation of the area with T is irregular 
for large values of T and the wing with A = 8 requires 
the larger values of T for spiral stability. 

The value of S f /s w required for. oscillatory and 

weathercock stability increases with a decrease in the 
wing chord (figs. 4 and 6). This effect is greater than 
the change in the actual fin area (fig. 6) and in p,, as 
shown in figure 2, and is primarily due to the decrease 
in the values of C- and C_ of the wing as compared 

with the fuselage size. It appears from these results 
that changes in fin area should be based on fuselage size 
and fuselage characteristics rather than on wing area. 
This conclusion is in agreement with the results from 
flight tests reported in reference 8. 
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11 feci of tail length on stability boundaries. - In- 
creasing the ratio of the tail length to the wing span 
i t /b decreases, almost directly, the value of Sf/S^ re- 
quired for stability (figs, 5 and 7). This decrease re- 
sults from the fact that the computed values of C n ^ 

vary directly with (Sf/S^) ^l t /b^ f but the closeness to 

a direct variation is also due -to the fact that the com- 
puted decrease in Cy^ is partly count er- balanced by the 

computed increase in 0 n . 

u r 

As the tail length is increased, the value of I\ 
and the possibility of obtaining stability, is increased. 
These effects are indicated by figure 7, where the effects 
of T are small but may be important for airplanes having 
large values- of k x /b, k 2 /b, and p,. 

GENERAL COMMENTS 



The aerodynamic characteristics of airplanes vary 
considerably; so choice of their values, their rates of 
change with various parameters, and their interference 
effects may greatly modify the magnitude of these results. 
The results of this investigation indicate general effects 
of changes in certain airplane parameters and do, not indi- 
cate the stability characteristics of a particular air- 
plane. 

The general results of this investigation may be 
summed up by pointing out the factors that make the at- 
tainment of lateral stability more critical and difficult: 

1. Large values of \x (high wing loading and high 
altitude) and large values of radii of gyration (weight 
distributed along the wings) 

2V Small wing chords 

3. Short tail lengths- (The magnitude of this ef- 
fect is comparatively small.) 

4. High lift coefficients and particularly the ne- 
cessity of obtaining stability both with and without flaps 
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Other factors,: such as the type and the degree of- stabil- 
ity desired and the choice of aerodynamic characteristics 
used, should have a large effect on the "interpretation of 
these results. 

Although the type and the degree of stability that 
affect the control and the riding qualities of the air- 
plane are outside the scope of this report, some of these 
factors should be considered. 

It is normally considered that, for satisfactory sta- 
bility characteristics, the oscillatory motion should be 
highly damped, that is, the airplane should have a large, 
amount of oscillatory stability. This stability can be 
obtained by proportioning the effective fin area and the 
effective dihedral angle so that their values will be 
positive, that is, when plotted on a diagram for a par- 
ticular airplane, the value will be in a stable area and 
well above and to the left of the zero oscillatory stabil- 
ity boundary of the airplane. The distance from the 
boundary, although an indication, is not a quantitative 
measure of the amount of damping. An airplane should be 
spirally stable, particularly for flying conditions of 
poor visibility. For good riding qualities in rough air, 
however, the amount of spiral stability should be small; 
spiral instability is generally considered to be prefera- 
ble to a poorly damped oscillatory motion. The value of 
c np should be positive, weathercock "stability; the 

amount, however, is dependent upon the spiral and the os- 
cillatory stability requirements. 



CONCLUSIONS 



Prom the analysis for lateral stability, assumed data 
for changes in certain airplane parameters being used, the 
following conclusions may be drawn: 

1. Increasing the values of the relative density of 
the airplane to the air, that is, increasing the wing 
loading or the altitude, makes the attainment of lateral 
stability increasingly difficult and critical, 

2. Increasing the radii of gyration more than cer- 
tain amounts makes the attainment of lateral stability 
very difficult and critical. 
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3. Increasing the aspect ratio of the wing by de- 
creasing the chord makes the attainment of lateral stabil- 
ity more difficult, 

4. Increasing the tail length, with the fin area in- 
versely proportional to the tail length (constant tail 
volume) , makes the attainment of stability slightly less 
difficult. 

5. The fin area required for lateral stability de- 
pends more upon the fuselage size than upon the wing area 
and, for a given fuselage, the required fin area should 
he increased slightly as the wing area is decreased. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , May 14, 1941. 
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Figure 1.- Variation of kj/b with kx/b fr o m design data available 

for pursuit- type airplanes and those used in the 
computations. 
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Fig. 3 
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Figure 3. Sffect of n on boundaries of spiral and oscillatory stability. A , 8 ; lt/b , 0.4 ; * x /b , 0.135; 
k z /b , 0.175 . 
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Fig. 3 
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12 16-4 O 4 

Effective dihedral angle, T, deg 
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Figure 3.- Kffeot of radii of gyration on boundaries of spiral and oscillatory stability. A , 8 ; X./b , 0.4 \ 
\i , 9.70 . * 
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Fig. 4 
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Effective dihedral angle, T, deg 
(a) C L , 0.3 (b) C L , 0.8 (c) C L , 1.4 (d) C L , 3.8 

Figure 4. Effect of wing chord on boundaries of spiral and oscillatory stability. l t /b , 0.4; k r /b , 0.135; 
k z , 0.175. 
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fig. 5a,b 
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figure 5a to d.- /Effect of tail length on boundaries of spiral and oscillatory stability. A , 8 ; k»/b , 0.125; 
kg/b , 0.175; |i ,9.70 . * 
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Figs. 6,7 
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